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ABSTRACT 

Tuo and thtaa dimanaional invitcid aolutioni for tha flou 
uxthin a ttanaonic aNial comptatior toXot at daiign ipaad ara 
comparad to laaat ananomatat maasuramanta av maximum tlou and 
naai stall oparating oointi Computational datails of tha 2-D 
axiiymmattic straam function solution and tha 3-D full Eular 
solution ara datcribad. Upitraam of tha rotor tha ?-D and 5-D 
solutions for radial distribution of ralativa Mach numbar and 
total prassura agraa uall uith tha data. Uithin tha bow uava 
systam and tna blada rou tha axisymmatric 2-D solution shous only 
gunlitativB agraamant uith tha data. Comparisons ara maria batuaan 
ralativa Mach mimbar contours, shock location, and shock 
strangth a maasurad and as ptadictad 'y tha 3-D coda. 
Comparisons at maximum flou shou taasonably good agraamant at all 
span-uisa locations but also soma disagraamant dua to viscous 
atfacts. Comparisons naar tha tip for tha naar stall casa shou 
axcallant agraamant Dounstrasm of the rotor tha inviscid 
computations agraa uith aach other but pradict hiahar prassura 
ratios than thosa maasurad. Eular codas typically raquira a 
dounstraam prassura as input. Sinc'< that prassura contro.ls tha 
computad mass flow and shock syst.-m. it aiust ba consistant uith 
an inviscid solution A proesduta for using an afficiant 
SMsymmatric coda to ganarata dounstraam prassura input for mora 
'.ostly Eular codes is discussed 


Computational and experimental 
techniques for the study of fl3us in 
turbomachinery are evolving concurrently 
and are mutually supportive. Experimet.ts 
provide data for empirical ccrrelations 
of phenomena that cannot yet be computed 
such as turbulence. Experimental data may 
also he used to verify direct 
comp..,tations of flou phenomena. 
Computational methods can cclve for 
thermodynamic quantities that cannot be 
measured by optical techniques. 
Computational methods are ilso becoming 
important tools for the development of 
advanced turh^machine components, and may 
be used to screen neu designs before 
resorting to more costly experiments. 

Advances in both computational and 
experimental fluid mechanics have been 
paced largely by advances in electronics. 
p.irticularly in digital data acquisition 
and processing. In a 1958 NACA report 
111. tip static pressures uere measured 
in a transonic compressor rotor using 
oscilliscope traces of signals from four 
pressure transducers embedded in the 
shroud. Expcrimantalists nou rely on 
digital data acquisition* and reduction 
systems coupled uith traversing pressure 
and temperature probes for pitchwise 
averaged flou measurements upstream and 
dounstream of rotors 121. 


Ti-'o-dimensional thro>*gh-flou analysis 
t,odes like IJI are the computational 
counterparts of these probe measurements, 
(’•jrrent data processing equipment coupled 
ijith electro-optical sensors allous the 
experimentalist to make detailed maps of 
the core flou in rotating compressors. 
19 and 5). Three-dimensional analysis 
codes like UI are the computational 
counterparts of current laser anemometer 
measurements. Experimentalists are gust 
beginning to take measurements inside the 

viscous layers on blade rous. Similarly, 
computational methods are nou being 
developed to analyse tuo-dimensior..il 
VISCOUS flous in cascades 17). It seers 
likely that the latest generation ot 
supercomputers and continuing advances in 
instrumentation uill allou both 
computational and experimental uork to 
progress even further, and that the 
resulting gains in understanding will 
lead to improvements in turbomachinery 
performance. 

It IS important to note that the 
early experimental and analytical uork is 
still extremely important. In fact this 
paper relies heavily on all of references 
(i through 6] mentioned above. Before 
going into the details of the present 
uork. It may be of interest to revicu 
current literature comparing optical 
measi.:rements and computations of flous in 
turbomachinery . 

Optical techniques that have been 
used for flou measurements in 
turbomachinery include holographic 
interferometry 181. the gas fluorescence 
technique 191. the laser-tuo-focus 
technique (L2F. also called laser transit 
anemometry or LTA) ( iO and 111. and laser 
anemometry (LA. also called laser Doppler 


vclocimctry or LDV) and 5|. Th« LA 
tachniqua raportad in 15] was used in tha 
prasant uork. 

Computational mathods usad to 
calculata flous in turbomachinery include 
2-D through- flou codaa IS# 12 > and 16 I ■ 

2- D steady blada-to-blada codas 

113 and 16), 2-D unsteady blada-to-blada 

codas (15 and 16), quasi-3-D codas 
(13 and 16), and fully 3-D codes 
(6 and 17). Tha 2-D through-flow coda in 
( 3 ) and tha 3-D coda in ( 6 ) ware usad in 
tha present uorlt. 

Tha following rafarancas compare 
optical maasuramants and computations of 
flous in turbomachinary. Several 
references have compared computations to 
L2F measurements in a transonic rotor 
designed and tested at DFVLR. Reference 
(12) compares these measurements with 
through-flow calculations at a full speed 
operating point. Reference (13) compares 
them with a quasi-3-D solution at a 
subsonic part speed operating point, and 
reference (16) compares them with a 
qvasi-3-D time-marching solution at three 
transonic operating conditions. A fully 

3- D so.' ution for transonic flow in the 

DFVLR rotor is presented in (17). 
Comparisons between LTA measurements and 
steady blade-to-blade solutions of the 
flow in a low speed four stage GE 

research rotor are given in (16). 
Reference (IS) gives an early comparison 
between LA measurements and a steady 
blade-to-blade solution in a GE transonic 
fan. Reference (18) compares gas 
fluorescence measurements with a 3-D 
time-marching solution for a HASA low 
aspect ratio transonic fan. 

The present work compares 
conventional probe measurements, LA 
irea«urements , 2-D through-flow 

computations, and 3-D Euler computations 
of transonic flous in a NASA-designed 
core compressor inlet rotor. These 
comparisons are made at design speed at 
two operating points, a maximum flow 
oint and a near stall point. Comparisons 
etucen measured and computed data are 
made for total pressures and relative 
tiach numbers along grid lines, for 
nominal relative hach numbers on the 
blade surfaces, for contour plots of 
relative Mach numbers and flow angles on 
grid planes, for shock locations and 
strengths, and for mass flous and total 
pressure ratios. These comparisons are 
intended to demonstrate the capabilities 
of the analysis codes for predicting two 
widely different flow conditions in an 
actual machine. 

In reference (8) similar comparisons 
were made between LA data and a 3-D Euler 
solution of the maximum flow operating 
point for this rotor. Good comparisons 
were found near the tip, but the computed 
shock structures at midspan and near the 
hub were incorrect. Bow waves were poorly 
resolved. Several improvements over 
reference (8) have been made in the 


resent uork. First, bow wave resolution 
as been improved by increasing the grid 
size and by clustering the grid around 
the leading edge. Second, shock 
structures have been improved by 
replacing measured exit hub static 
pressure input to the code with values 
calculated using a 2-D through-flow code. 
The procedure used to calculate 
appropriate boundary values is described 
in detail later. Third, the near stall 
results presented here may be of more 
interest to designers since these results 
are . more representative of a typical 
compressor design point. 


EXeSRIt lg NT M r APPARATUS 
Compressor Rotor 

Rotor 33 was designed at NASA Lewis 
Research Center as an inlet rotor for a 
core compressor. The floupath consists of 
a constant radius hub and a converging 
tip. At the inlet the tip radius is 
256 mm (10 in) and the hub/tip radius 
ratio is 0.7 The rotor has 52 blades 
with a tip chord of 66.5 mm (1.75 in) and 
a tic solidity of 1.68 .At the design 
speed of 16,100 rpm the tip speed is 
626 m/sec (1398 ft/sec.) 

The rotor was tested without inlet 
guide vanes or stators, thereby 
eliminating any circumferential 
variations that would be introduced by 
stationary blade rows. Radial surveys of 
total pressure, total temperature, flow 
angle, and static pressure were taken 
upstream and downstream of the blade row 
using conventional probes (2). The 
upstream and downstream measurement 
stations, designated stations one and 
two, were located one chord upstream and 
one-half chord downstream of the rotor 
respectively . 

Data tor this paper were taken at 
design speed at two throttle settings, a 
wide open maximum flow setting and a 
partially closed near stall setting. At 
maximum flow the rotor pressure ratio was 
1.69 and the mass flow was 208 kg/sec/m^ 
(62.60 Ibm/sec/f t> ) . Mass flows were 
measured using a calibrated orifice 
located far upstream. Inlet relative 
Mach numbers at maximum flow ranged from 
1.16 near the hub to 1.35 near the tip. 
Near stall the rotor pressure ratio was 
1.66 and the mass flow was 191 kg/sec/m^ 
(39.12 Ibm/sec/f t* ) . Inlet relative Mach 
numbers near stall ranged from 1.05 n«>ar 
the hub to 1.32 near the tip. 


Laser Anemometer 

The fringe type laser anemometer 
(LA) ( 6 ) IS a single channel dual beam 
system with on-axis backscatter light 
collection. A dedicated minicomputer, a 
disk storage system, and a CRT display 
were used for real time data acquisition 
and processing. The minicomputer was also 
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connacted to a larga central computer for 
additional data processing and graphical 
output. 

k window fabricated from commercial 
window glass provided optical access to 
the compressor. The window conformed to 
the outer endwall contour and measured 
102 mm axially by 51 mm in circumferenco 
(11 degrtes arc) by 1 mm thick. 

Seed particles consisted of 
spray-atomized rhodamine 6G dye dissolved 
in a solution of benzyl alcohol and 
ethelyene glycol. The particles were 
injected into the flow A60 mm upstream of 
the rotor. Exposed to laser light the 
seed particles fluoresce orange. An 
orange-pass optical filter placed in 
front of the photomultiplier selectively 
removed unwanted green laser light 
reflected from the hub and blades. The 
seed particle diameter was determined by 
equating the particle velocity lag 
measured downstream of the rotor passage 
shock to tnat predicted oy a 
one-dimensional Stokes drag model and 
solving for the unknown particle 
diameter. The particle diameter found by 
this procedure was 1.2- 1.5 micrometers 
(51. 

An efficient data acquisition system 
made it possible to map the 
blade-to-blade distribution of a velocity 
component at a given axial and radial 
position very quickly. Whenever a seed 
particle crossed the probe volume the 
particle velocity and the rotor shaft 
angular position ware recorded as a data 
pair. Approximately 50 measurements were 
made at each of 1000 different shaft 
ositions distributed as 50 positions per 
lade passage over 20 consecutive blade 
passages. The 50 velocities and flow 
angles at each shaft position were 
averaged to give time-averaged velocity 
distributions over 20 consecutive 
passages. The time-averaged distributions 
were analyzed to detect variations 
between individual passages (5). The 20 
time-averaged velocity distributions were 
spatially-averaged to yield the velocity 


distribution across an average passage. A 
typical run consisting cx 50.000 
measurements at a single axial and radial 
position took between 15 and A5 seconds. 

Errors in the LA measurements arise 
from a number cf sources. The sources of 
error and the corresponding error 
magnitudes have been discussed in detail 
in [ 5 I and will only be summarized here. 
Passage-to-passage velocity variations 
were on the order of five percent for 
most runs. Spatial averaging of data 
obtained across 20 blade passages 
improves the statistical accuracy of the 
data. Compressor speed drift during a run 
was on the order of 0.5 percent and the 
repeatability of a given run was one 
percent. The error in velocity and flow 
angle measurements which encompassed 95 
percent of the data was generally five 
percent or less except in regions 
immediately downstream of shocks where 
seed particle velocity lag was the 
dominant error. The observed lag distance 
required for the seed particles to 
achieve 95 percent of the post-shock gas 
velocity was found to be 15 percent of 
chord for the maximum flow case at 15 
ercent span from the tip. The LA data 
as not been corrected for particle lag 
since velocity variations due to particle 
lag cannot in general be distinguished 
from those occurring due to the 
compression process. 

Figure 1 shows tl.e LA measurement 
stations and probe survey stations 
superimposed on the finite difference 
grid used with the 5-D code. Typically 
22 measurement stations lie along each of 
three design streamlines located at 15. 
50. and S5 percent span from the tip with 
12 of the axial stations within the blade 
row. The three radial locations lie 
within one percent span of the grid lines 
shown. At the near stall operating point 
LA measurements were only taken at 15 
percent span from the tip. 


COMPUTATIONAL METHODS 

Tuo-Dimensional Through-Flow Code MERIDL 

The through-flow analysis code used 
in the present work is the MERIDL code 
(3 1. The code solves an axisymmetric 
inviscid momentum equation for a stream 
function along a midchannel hub-shroud 
stream surface. The stream function is 
defined to satisfy continuity 
identically, and is differentiated 
lumerically to find individual velocity 
components. Thermodynamic quantities are 
found from isentropic relations. 
Rotation, passage convergence, blade row 
turning, and blockage are all accounted 
for it- the formulation. A successive over 
relaxation finite difference technique is 
used to solve the nonlinear stream 
function equation. An orthogonal finite 
difference grid with 10 points radially 
and AO points axially (10 points upstream 
and downstream, 20 within the blade r^jw) 
was used for the present calcuJ . ■ ,us . 
While the grid is not shown here, had 
the same extent as the 3-D grid shown in 
Figure 1. Input to MERiDL consists of 
blade, hub, and tip geometry, mass flow, 
upstream total conditions, up and 
downstream whirl (rVe) distributions, and 
axial and radial distributions of total 
pressure loss. Upstream total conditions 
were input as constants radially, and the 
upstream whirl was set to zero. The exit 
whirl was taken as a smooth curve fit 
through the measured probe data, 
excluding points obviously in the endwall 
boundary layers. 

Three-Dimensional Euler Code BLADE3D 

The 3-D analysis code used in the 
present work is the BLADE3D code 
described in 16). Only a brief 
description is given here. The equations 
solved are the inviscid unsteady Euler 
equations written in cylindrical 
coordinates and in non-conservative form. 
The equations include the continuity 
equation, the axial, radial, and 
circumferential momentum equations, and 
the energy equation, which are solved for 
the density, three velocity components, 
and total energy. Pressure is found from 
the ideal gas law. Blade row rotation is 
: peci f led . 

MacCormack's explicit time-marching 
finite difference method in split 
operator •‘orm (19) is used to advance the 
unsteady solution in time from an initial 
guess to a converged steady solution. The 
metliod IS second order accurate and 
conditionally stable. Shocks are captured 
automatically and smeared over several 
grid points. Artificial viscosity terms 
are added to improve stability near 
shocks, but they also increase shock 
smearing . 

Algebraic mappings are used to 
transform the complex flowpath geometry 
into a rectangular computational domain. 


The computational mesh had 100 axial 
points (32 points upstream and 
downstream, 3A within the blade,) 17 
circumferential points, and 18 radial 
points. Axially the grid points were 
clustered around the leading and trailing 
edges to impove resolution of bow and 
shock waves. Points were stretched 
upstream and downstream to allow 
imposition of axisymmetric far field 
conditions. Radially and 
circumferentially the points had constant 
spacing. A sheared blade-to-blade grid 
was used. Figure 1 shows a meridional 
view of the grid. 

At the inlet grid plane the flow is 
assumed to he axisymmetric and is matched 
to desirud conditions at upstream 
infinity where total conditions are 
specified, the radial velocity is zero, 
and the whirl is specified (here also 
zero.) The non-ref lective inlet boundary 
condition is based on the one-dimensional 
method of characteristics. The derired 
axial velocity and static temperatur * at 
upstream infinity are combined to fo-m a 
single variable, the value of the 
downstream-running characteristic. This 
value may be computed from the desired 
mass flow or upstream Mach number using 
isentropic rela<-*ons, or from MERIDL 
output. Although the BLADE3D solution is 
sensitive to the input characteristic 
value, all methods of computing that 
value give consistent results. The 
intersection of the specified 
downstream-running characteristic and the 
computed upstream-running characteristic 
at the inlet determine the inlet static 
temperature and axial velocity. Thus the 
mass fElow is computed as part of the 
solution and may not necessarily match 
the measured value. 

At the exit grid plane the flow is 
also assumed to be axisymmetric and the 
radial velocity is set to zero. The exit 
circumferential velocity, density, and 
downstream-running characteristic are 
extrapolated from within the flow field. 
Only one boundary value is specified, the 
hub static pressure at the exit of the 
computational domain. A radial 
equilibrium equation is integrated from 
hub to tip at the exit to determine the 
radial pressure distribution, which 
supplies the final unknown at the exit. 

Other boundary conditions used in 
the BLADE3D code include periodicity up 
and downstream of the blade row, tangency 
on the hub, shroud, and blade surfaces, 
and a Kutta condition specifying constant 
static pressure across the blade trailing 
edge . 

Exit Pressure Specification 

Euler codes typically require that 
the exit static pressure be specified as 
a boundary condition when the exit flow 
IS subsonic. This input variable 
controls the computed mass flow and siiock 



structura, and must ba specified 
correctly to produce reasonable 
solutions. Experienca with the BLADE3D 
coda has shoun that use of measured 
values of exit hub static pressures tends 
to produce solutions with hig.iur mass 
flous and stronger shocks than t.iose 
measured. This appears t» be due to the 
neglect of viscous blockage and losses. 
Without these viscous effects, the 
BLADE3D code predicts higher mass flous 
for a oiven static pressure rise, or 
conversely, predicts a higher static 
pressure rise for a given mass flou than 
uould be found in a viscous flou. Hence, 
to compute a cose uith a given mass flou 
a designer must input an exit static 
pressure someuhat higher than measured. 
That exit pressure can be varied 
Iteratively to match the desired mass 
flou, but computational times may become 
prohibitive. 

Since the desired mass flou is an 
input variable for tha MERIDL code and 
static pressures are included in the 
output, and since a MERIDL solution can 
be run in about a minute of computer 
time, the MERIDL code can be used 
effectively to estimate exit hub static 
pressures for the BLADE3D code. 

Measured values of mass flou, 
upstream total conditions, and inlet and 
exit uhirl are input to MERIDL as 
desf.tibed earlier. One piece of input 
remains to be described, the axial and 
radial distributions of total pressure 
loss. Since MERIDL is incapable of 
predicting either shock or viscous 
losses, tne amount of loss specified can 
make a MERIDL solution compare uell uith 
either measured data or a BLADE3D 
solution. For best comparison uith 
measured data, measured losses must be 
input. But for best comparison uith a 
BLADE3D solution, only shock losses must 
be input. 

The axial distribution of the loss 
has little or no effect on the exit 
static pressure, so losses may be 
distributed evenly uithin ^h« blade rou. 
Radial distributions of shock losses in 
rotor 33 uere calculatsd as part ox the 
probe data reduction u:ing the 
Miller-Hartmann model described in (1). 
The model estimates shock losses to be 
those occurring across a single normal 
shock standing near tha passage entrance. 
The Mach number ahead of the shock is 
taken a? the average of the inlet 
relative Mach number and a someuhat 
higher value on the suction surface near 
the shock intersection. 

At the near stall operating point 
rotor 33 has a single shock system like 
the one in the Miller-Hartmann model, but 
at the maximum flou point it has a tuo 
shock system that is not modalled 
correctly. Overall the BLADE3D near stall 
solution compares better uith the LA data 
than the maximum flou solution does. Much 
of this IS attributed to tha better fit 



a(kf /tac/B^) 

aClbe/erc/f t^) 

tot .praa. ratio 

u 

206. 

42,6 

1.49 

BLADE 3D 

2iB. 

44.6 

1.71 

SRIDL 

208.* 

42.6 * 

1.64 


Table l. Comparison of measured and 
computed mass flous and total pressure 
ratios, maximum flou point. 


**Mass flou is set explicitly as an input 
to MERIDL. 


of the near stall shock structure to the 
Miller-Hartmann model. 

Computational Details 

Initial conditions for the maximum 
flou case assumed constant upstream 
velocity based on a specified inlet Mach 
number, velocity variations uithin the 
blade rou based on blade turning, and 
constant velocity dounstream. Solutions 
uere run on an IBM 3030 AP computer. The 
time step used ranged from 7S to 90 
percent of the stability limit. 
Convergence to a steady state uas 
determined by monitoring the trailing 
edge pressure distribution uhich reflects 
the Kutta condition and the blade 
loading. Trailing edge pressures uere 
converged to the fourth significant 
figure . 

The maximum flou case took 
approximately 12 CPU hours (A500 time 
steps) to converge. The near stall case 
uas calculated as the time-accurate 
response of the maximum flou case to a 
9.6 percent increase in exit hub static 
pressure and a change in inlet conditions 
uhich resulted in an inlet Mach number 
drop of 0.06. This case took nearly 
tuice as long to converge as the maximum 
flou case did because of the additional 
time necessary for the dounstream 
pressure rise to travel upstream against 
a high sub.sonic flou. Starting from the 
nearly constant initial conditions 
described earlier, the near stall case 
converges about as fast as the maximum 
flou ca.se. Thus, uhen only a steady 
solution IS desired, it is often 
computationally more efficient to start 
uith a simple initial guess than uith a 
converged solution for a different flou. 

The computer times quoted above are 
large. It should be noted that they uere 
accumulated over many runs by restarting 
the code. Computer times can be reduced 
substantially by reducing grid size and 
grid clustering. 


umis m piscussiQH 

Maximum Flou Operating Point 

Since mass flou is computed by 
BLADE3D and pressure ratio is computed by 
both BLADE3D and MERIDL. these global 
arameters are of particular interest, 
able 1 shous that the mass flou 
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calculated by BLADE3D is ^.8 percent 
higher than the measured orifice value. 
This IS a substantial error at this 
operating point. Nevertheless. other 
aspects of the BLADE3D solution agree 
reasonably uell uith the LA data. Total 
pressure ratios calculated by both codes 
are high; BLADE3D is 14.8 percent high 
and nERIDL is 10.0 percent high. 
Differences betueen the BLADE3D and 
MERIDL solutions are attributed to 
differences in shock losses calculated by 
BLADE3D and by the tliller-Hartmann model. 

Radial distributions of 
circumferentiallv-averaged total pressure 
and relative Macn number arc compared to 
probe measurements upstream and 
dounstream of the rotor in Figures 2 and 
3. Axial locations labeled station I 
and station 2 are shotin in Figure 1 . 
Upstream (station 1) the calculations and 
measurements arc in good agreement except 
for the error due to the high BLADE3D 
mass flou. Dounstream (station 2) the 
high total pressures calculated 
inviscidly are evident. Enduall boundary 
layer effects that cannot be predicted 
inviscidly are evident in the probe 
measurements, especially at the tip of 
station 2. 
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figure 4. • Alibi disUlbulion of relitiv* Mbch num- 
btr, IS ptrcenf totn from the tip, mM-qtp. rrviii- 
mum flow point. 


Axial distributions of relative Tach 
number at 15 percent span from the tip 
are compared in Figure 4 The flERIDL 
calculations are of course axisymmetric . 
The BLADE30 and LA results are shoun 
along a mid-gap line. The flERIDL results 
shou little upstream influence of the 
blades, uhile the BLADE3D and LA results 
shou relative Mach number fluctuations 
due to a bou uave system. The computed 
bou uave is weaker than that measured. 
Uithin the blade rou the measured and 
BLADE3D computed flous accelerate 
slightly, decelerate across an oblique 
passage shock, and finally decelerate 
sharply across a normal shock near the 
trailing edge. The calculiited normal shock 
IS stronger than the .measured shock, 
possibly due to the error in the computed 
mass flou or to the neglect of viscous 
blockage. Discrepancies betueen the 
shock }umps shoun and the normal shock 
tables are due to the inclination of the 
rear shock relative to the line along 
uhich the results are shoun. The flERIDL 
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results shou surprising agreement inside 
the blade rou. but the agreement is 
coincidental. The initial acceleration in 
the HERIDL results is due to blade 
blockage. The rapid deceleration that 
appears to be a shock is due to rapid 
turning of the flow to match the 
specified dounstream uitirl. Dounstream 
of the blade rou both computed flous 
reaccelerate slightly due co tip 
convergence. 

Circumferential distributions of 
relative nach number at IS percent span 
from the tip and at 78 percent chord 
(through the rear shock) are shoun in 
Figure 5. The abscissa ox the plot 
shouing percent gap is measured from the 
suction surface of a leading blade to the 
suction surface of the follouing blade. 
At this axial location the flou near the 
.'(uction surface (left) is upstream of the 
rear shock but the flou near the pressure 
surface (right) has prssed through the 
shock. It IS possible that the measured 
drop off in Tlach number near the pressure 
surface is caused by shock-boundary layer 
interaction. The shock location is 
clearly defined in both the LA 
measurements and the BLADE3D calculations 
as an abrupt drop in Hach number near 
mid-gap. Shock smearing evident in the LA 
measurements is du* to seed particle lag; 
and shock sm(»aring in the BLADE5D 
calculations is due to artificial 
viscosity. Blade boundary layer effects 
may be evident in the measurements. 
riERlDL results look reasonable as an 
axisymmetric average. MERIDL results are 
omitted from the remaining maximum flou 
results. 
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Tabic 2. Compariton of measured and 
calculated shock loss parameters, rear 
shock, mid-gap, maNimum tlou point. 


Blade surface relative Mach number 
distributioi's along the design streamline 
at 15 percent span from the tip are shoun 
}n Figure 6. LA measurements of surface 
relative Mnch number uere taken from 
plots like Figure 5 as the first 
circumferential point for which a 
statistically significant number of 
measurements was taken (at least 20 
measurements ) At a tip chord Reynolds 
number of about 6*10*. an unseparated 
turbulent blade boundary layer thickness 
would be of the order of the measurement 
point spacing; so it is felt that the 
first statistically significant point 
should be near the edge of the core flow. 
Nevertheless, viscous effects may be 
present in the data. The calculations 
and measurements show only a qualitative 
agreement in shock location and the 
rearward blade loading. 

Figure 7 shows contours of computed 
relative hach number on a mean flow 
surface at mid-chsnnel. and is included 
to show the three-dimensional nature of 
the solution. Most of the upstream and 
downstream regions have been omitted. 
Dashed lines indicate shocks. The rear 
passage shock is slightly curved from hub 
to tip. Near the hub it is difficult to 
distinguish whether there is a rear shock 
or 3 ust a diffuse compression. 

Computed and measured blade-to blade 
relative Mach number contours are 
compared side by side at the three 
span-wise locations in Figure 8. Shocks 
are shoun by dashed lines. Computed shock 
locations were determined from axial 
plots of relative Mach number like Figure 
•t as grid points with maximum Hach 
numbers before rapid drops. These jrid 
points were located on the blade-to blade 
plots and fit with a smooth curve. 
Experimental shock locations were 
determined similarly except that 
Circumferential pTots like Figure 5 were 
used. At 85 percent span from the tip the 
computations and measurements are in 
excellent agreement. Both show a bow wave 
with an exit Mach number of I . 2 The 
calculations show a stronger suction 
surface acceleration than the data. Both 
also show a forward passage shock with 
similar exit Mach numbers and small sonic 
regions on the pressure surface near the 
leading edge. Finally both show a diffuse 
compression to subsonic with no obvious 


rear shock. The calculated flow exits at 
a considerably lower Hach number than the 
measured flow. At midspan the comparison 
IS similar except that here both show a 
definite rear shock. The measured rear 
shock IS forward in the passage relative 
to the calculated shock, and may even be 
a reflection of the front shock. At 15 
percent span from the tip the LA 
measurenents show a stronger bow wave and 
passage shock than the BLADC3D solution . 
The calculations and measurements both 
show very strong rear shocks attached to 
the trailing edge. These results show 
improved bow wave resolution and rear 
shock location over the results in 
reference (B| due to grid refinement and 
adjustment of the exit hub static 
pressure . 

Shock angles and total pressure loss 
coefficients for the forward passage 
shock at the three span-wise loc-cions 
and at mid-gap are tabulated in Table 2. 
The ’shock inclination angle. <*., is 
measured from the axial direction. The 
shock loss coefficient is defined in [1) 
as : 

uTi»(l-P, /Pj )/(l-p, /Pi ) 

where upper and lower case letter p refer 
to total and static pressure 
respectively. The total pressure ratio 
P, /Pi is obtained from the normal shock 
relations using the component of relative 
Mach number ahead of and normal to the 
shock. There is some disagreement between 
the LA and BLADE5D shock loss 
coefficients, especially at the hub. Part 
of the disagreement is due to small 
d.vscrepancies ahead of the shock; part is 
due to shock inclination in the 
meridional plane; but most is due to 
difficul les in determining the shock 
angle from stiock locations found as 
described earlier. 




■ (kg/kcc/a^) 

ii( Iba/scc/f t 

tot .pret. ratio 

.JL 

191. 

J9.1 

1.66 

BLADE3D 

188. 

38.5 

1.866 

HERIDL 

191. * 

39.1* 

1.864 


Table 3. Comparison of measured and 
calculated mass flous and total pressure 
ratios, near stall point. 


*f1ass flow is set explicitly as an input 
to MEItlDL. 
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.“liar Stall Operating Point 

nass flow and total pressure ratio 
fo: the near stall case are compared in 
Table 3. The BLADE3D mass flow is 1.5 
percent low . and the BLADEJD and MERIDL 
total pressure ratios are both 12. A 
percent high, probably due to neglect of 
viscous losses and blockage. 

Hadial distributions of 
circumferentially-averaged total pressure 
and relative Mach number shown in Figures 
t and 10 arc similar to the maximum flow 
case except that here the BLAOEJD and 
nEi'IOL results are in better agreement. 


Axial distributions of relative Mach 
number at 15 percent span from the tip 
are shown in Figure 11. The BLADE3D Mach 
numbers show good agreement with the LA 
measurements in the strong bow wave 
system and in the location of the passage 
shock. The passage shock jump is over 
predicted, again possibly due to the 
neglect of viscous blockage. The MEPIDL 
results behave as before, shoiang an 
initial acceleration due to blockage 
followed by a rapid compression due to 
turning. MEPIDL cannot capture the 
upstream wave system or the passage 
shock, yet the MEPIDL results agree 
fairly well outside the blade row. 
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Circuihferential distributions of 
relative Mach number at 15 percent span 
from the tip and at 25 percent chord 
(through the passage shock) are shoun in 
Figure 12. At this oKial location the 
flou near the suction surface (left) is 
urstream of the passage shuck but the 
flou near the pressure surface (right) 
has passed through the shock. The Mach 
number distributions shou slight 
discrepancies in shock location and 
strength. nCRIDL results agree uith flou 
condition', on the upstream side of the 
shock. M'JRIDL resu.lts are omitted from 
the remain.' nq near stall results. 

Fioure 15 compares computed and 
measured blade surface relative Mach 
numbers at 15 percent span from the tip. 
Suction surface shock locations agree to 
uithin the measurement spacing. The 
computed near stall blade loading is 
dramatically shifted to the front of the 
blade compared to the maNimum flou 
solution in Figure 6. 

Computed relative nach number 
cont irs on a mid'Channel mt^an-flou 
surface in Figure 1^ shou a complete 
change in shock structure from the 
maKimum flou solution in Figure 7. In 
the near stall solution the Sou uave and 
passage shock are stronger and straighter 
than in the maximum flou solution. The 
rear passage shock has been eliminated. 



figurt 15. - Ama fltstnbulion 01 rtlitiv* M*ch nu("t«r on 
tht t;«d» lurtKtl, 15 pt( c»nt spn lro(" m» lip, ntir 
st«ll potr.' 
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T«bl» ^ . Comparison of measured and 
calculated shocK loss parameters > passage 
shock, mid-gap, near stall point. 


The shock structure at 15 percent 
span from the tip is shown in Figures 15 
and 16 which compare computed and 
measured contours of relative Mach 
numbers and flow angles. Unlike the tip 
Mach number contours for the maKimum flow 
case in Figure 8, Figure 15 shows e 
pronounced bcu wave and passage shock 
system uiv:h excellent agreement between 
the BLADE3D and LA results. Computed 
relative Mach numbers leaving the shock 
are low, however. The wake-like contours 
leaving the trailing edges in the 
computed results show an inviscid 
trailing edge slip line. Relative flow 
angle contours in Figure 16 shew an 
upstream angle of 68 to 69 degrees. The 
suction surface metal angle at the 
leading edge is about 6<i degrees. The 
high incidence incoming flow turns 
smoothly through expansion fans 
originating one or two blades away to 
become well aligned with the suction 
surface at tht leading edge. 

Table 9 compares the measured and 
computed shock loss parameters. Since the 
passage shock is well defi.ned in this 
case, the shock location and angle are 
relatively easy to determine. Slight 
discrepancies in flow conditions ahead of 
the shock account for a small 
disagreement in the shock loss 
parameters . 


CONCLUDING REMARKS 

Current laser anemometer (LA) 
measurements are of sufficient accuracy 
for verification of inviscid codes. Shock 
smearing in the LA measurements due to 
seed particle lag is of roughly the same 
magnitude as numerical shock smearing, 
and both .leed to be reduced. Either a 
better method of estimating boundary 
layer edge velocities or a means of 
directly measuring blade surface 
pressures is needed. More detailed 
measurements of blade and endwall 
boundary layers will be needed for 
verification of viscous codes. 

Inviscid calculations from an 
axisymmetric through-flow code (MERIDL) 
and a three-dimensional Euler code 
(BLADE3D) compared well with probe data 
upstream of the blade row. Downstream, 
total pressures are high and relative 
Mach numbers are low due to the neglect 
of viscous effects. Within the blade row 
the axisymmetric solutions agree only 
qualitatively with the other results. 


The BLA0E3D code accurately 
predicted Mach number and flow angle 
distributions and shock structure within 
an axial compressor rotor at maximum flow 
and near stall operating points. The near 
stall solution agrees better with LA 
measurements than the maximum flow 
solution does. This is probably because 
most of the flow turning in the near 
stall case occurs across a forward 
passage shock away from viscous effects . 
But in the maximum flow case a second 
normal shock near the blade trailing edge 
(where viscous effects would be greatest) 
controls the blade row pressure rise. 
Present results for the maximum flow case 
have improved on the earlier results 
given in reference (8J. Bow wave 
resolution was improved through grid 
clustering and rear shock structure was 
improved by adjusting the exit hub static 
pressure. 

Downstream pressure boundary 
conditions typically needed by Euler 
codes are difficult to specify. Measured 
downstream pressures produce erroneous 
mass flows and shock structures because 
measured pressures contain viscous losses 
not accounted for in inviscid codes. 
Computational times may prohibit 
iterating on boundary conditions. An 
axisymmetric through-flow code like 
MERIDL may be used as a consistent and 
computationally efficient tool for 
choosing downstream boundary values. 
MERIDL input includes upstream total 
conditions, mass How, and whirl, and 
downstream whirl, all input as measured. 
Shock losses must be specified and may be 
estimated using the Miller-Hartmann model 
along with measured probe data. Loss 
distribution through the blade row does 
not significantly effect the computed up 
or downstream solution. Downstream 
pressures calculated in this manner are 
higher than measured pressures but are 
appropriate boundary values for Euler 
cudes, producing the best agreement 
between mass flows and shock structures 
measured in a viscous flow and computed 
inviscidly. In future three-dimensional 
viscous codes the measured downstream 
pressure should suffice. 
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